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Abstract 

This  document  constitutes  a  final  report  for  our  research  in  electromagnetic  scatter¬ 
ing  during  1997.  The  technical  research  accomplishments  divide  conveniently  into  the 
areas  of  time-domain  electromagnetics,  high-frequency  diffraction,  and  surface  integral 
equations.  As  principal  investigator.  Dr.  Arthur  D.  Yaghjian  carried  out  the  major 
portion  of  the  research,  and  coordinated  his  efforts  with  those  of  T.B.  Hansen  and  R.A. 
Shore  of  Rome  Laboratory,  R.A.  Albanese  of  Brooks  AFB,  S.W.  Lee  of  DEMACO,  and 
R.  Tolimieri,  M.  An  and  A.J.  Devaney  of  A.J.  Devaney  Associates. 


1  Technical  Summary  of  Accomplishments 

In  1997  our  research  accomplishments  in  electromagnetic  scattering  divide  conveniently  into 
three  main  areas: 

•  Time- Domain  Electromagnetics 

•  High-Frequency  Diffraction 

•  Surface  Integral  Equations. 

1.1  Time- Domain  Electromagnetics 


Plane- Wave  Theory  of  Time-Domain  Fields 

In  the  area  of  time-domain  electromagnetics.  Dr.  Arthur  D.  Yaghjian  and  Dr.  Thorkild  B. 
Hansen  have  been  working  together  to  develop  the  “plane-wave  theory  of  time-domain  fields,” 
and  to  apply  this  theory  to  the  near-field  scanning  of  antennas  and  scatterers.  Previously, 
near-field  measurement  techniques  had  been  formulated  only  for  antennas  and  scatterers 
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operating  over  a  narrow  frequency  band,  that  is,  only  in  the  frequency  domain.  The  new 
time-domain  formulation  allows  powerful  near-field  measurement  techniques  to  be  extended 
to  pulsed  radiators  (antennas  or  scatterers)  operating  over  arbitrarily  wide  bandwidths.  This 
time-domam  near-field  measurement  formulation  was  published  in  the  IEEE  Transactions 
on  Antennas  and  Propagation  [1]  and  received  the  1995  Schelkunoff  Prize  Paper  Award. 

During  1997  Yaghjian  and  Hansen  expanded  and  refined  their  “plane-wave  theory  of 
time-domain  fields,  and  completed  about  90%  of  a  book  under  this  title  [2]  scheduled 
to  be  published  by  the  lEEE/Oxford  University  Press  in  1998.  Numerous  additions  and 
improvements  were  made  to  the  theory.  For  example,  Yaghjian  found  a  direct  method 
of  deriving  the  dispersion  relations  for  general  electromagnetic  materials  that  reveals  the 
natural  sufficient  conditions  on  the  constitutive  parameters  required  for  the  validity  of  the 
dispersion  relations.  This  derivation  is  given  here  in  Appendix  A. 

o«  of  time-domain  far  fields  was  submitted  to  Amerem- 

96,  and  this  year  an  expanded  version  of  the  paper  was  published  as  a  book  chapter  [31. 
Recently,  Yaghjian  has  been  invited  by  the  special-issues  editor  of  the  Proceedings  of  the 
IEEE  to  submit  a  review  article  on  time-domain  near-field  scanning.  He  has  also  been  invited 
to  subinit  a  paper  on  ultra-wideband,  short-pulse  electromagnetics  to  Euroem-98,  to  be  held 


During  the  year  we  have  assisted  Analytic  Designs  (Columbus,  OH)  in  their  implemen¬ 
tation  of  planar  near-field  measurements  in  the  time  domain.  We  have  also  been  a^ked  by 
Satimo,  Inc.  (Acworth,  GA)  and  Supelect  (Paris,  France)  to  work  with  them  in  their  com- 

a  time-domain  spherical  near-field  measurement  system.  In  addition, 

Near  MM  ^BIR  evaluation  committee  at  Rome  Laboratory  chose  “Time-Domain 
Near-Field  Measurements”  as  their  top  priority  Phase  II  SBIR  topic 

1  i"  Appendix  B.  Final 

Iff.  d  hapten  and  sect.ons  2.3-2.4  of  the  book  have  been  completed.  As  part  of 
the  effort  to  develop  and  document  the  plane- wave  theory  of  time-domain  fields.  Yaghjian 
ansen,  aong  with  Professor  A.J,  Devacey,  discovered  and  published  as  an  lEEE-APS 
Letter  a  profoundly  simple  method,  using  plane-wave  spectra,  for  determining  the  position 

S  t  t  Transactions  on  Antennas  and  Propagation  has  asked  us  to  submit 

whi!h  U  o Iffl  d  T  ?  application  of  this  method, 

il  V  v«  f  li  "!i  "ictccmining  the  minimum  source  region,  and  the  limits 

g  es  for  the  smallest  possible  sizes  of  antennas  that  can  radiate  specified  far-field  patterns 
provide  antenna  engineers  with  a  new  set  of  practical  design  tools. 


Trajectories  and  Fields  of  Bound  and  Unbound  Charges 

The  1992  book  [5]  by  Yaghjian  has  received  notable  acclaim  by  the  physics  community 
For  example,  Pro  essor  Ementus,  Fritz  Rohrlich  of  Syracuse  University,  considered  by  man^ 
physicists  as  the  foremost  expert  on  the  classical  equations  of  motion  for  charged  particles^ 
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writes 


This  IS  a  remarkable  book.  Arthur  Yaghjian  is  by  training  and  profession  an 
electrical  engineer;  but  he  has  a  deep  interest  in  fundamental  questions  usually 
reserved  for  physicists.  Working  largely  in  isolation  he  has  studied  the  relevant 
papers  of  an  enormous  literature  accumulated  over  a  century.  The  result  is  a 
fresh  and  novel  approach  to  old  problems  and  to  their  solution. 

Physicists  since  Lorentz  have  looked  at  the  problem  of  the  equations  of  motion 
0  a  charged  object  primarily  as-a  problem  for  the  description  of  a  fundamental 
partmle  typically  an  electron.  Yaghjian  considers  a  macroscopic  object,  a  spher¬ 
ical  insulator  with  a  surface  charge.  He  was  therefore  not  tempted  to  take  the 
point  limit  and  he  thus  avoided  the  pitfalls  that  have  misguided  research  in  this 
held  since  Dirac’s  famous  paper  of  1938. 


Perhaps  the  author’s  greatest  achievement  was  the  discovery  that  one  does  not 
need  to  invoke  quantum  mechanics  and  the  correspondence  principle  in  order 
o  exclude  the  unphysical  solutions  (runaway  and  pre-acceleration  solutions). 
Rather  he  discovered,  the  derivation  of  the  classical  equations  of  motion 
from  the  Maxwell- Lorentz  equations  is  invalid  when  the  time  rate  of  change  of 
-he  dynamical  variables  is  too  large  (even  in  the  relativistic  case).  Therefore, 
solutions  that  show  such  behavior  are  inconsistent  consequences.  The  classical 
theory  is  thus  shown  to  be  physically  consistent  by  itself.  It  is  embarrassing-to 
say  the  least— that  this  observation  had  not  been  made  before. 


This  work  is  an  apt  tribute  to  the  centennial  of  Lorentz’s 
in  which  he  first  proposed  the  Lorentz  force  equation.” 


seminal  paper  of  1892 


RohHfrbl  rf  i  K  terminology  of 

fcce  erat  d  b”''^  ^hen  an  electron  is 

accelerated  by  an  incoming  electromagnetic  pulse,  the  electron  experiences  a  “radiation  re- 

the  ve  externally  applied  force.  In  1997  we  began  work  looking  into 

e  version  of  the  Lorentz-Abraham-Dirac  equation  of  motion  derived  by  Yaghjian  [5  eq 
8.76]  for  an  accelerated  charged  particle  ^  ^  ^ 


pi  2  ^  /V 


ci  u  .duj  du- 


t':  corresponding 


^ext(^)  —  TTl 


d(7u) 


bireoc^ 


+  ^2  |up  +  ^(u-u)2  uj  (2) 


where  u,  e  and  m  are  the  velocity,  charge  and  mass  of  the  particle.  Of  course,  t  is  the  time, 
c  is  the  speed  of  light,  cq  is  the  permittivity  of  free  space,  Fej;t(t)  is  the  externally  applied 
force,  and  7  =  (1  —  u^/c^)~2.  The  force  is  applied  at  t  =  0  and  T]{t)  is  a  function  that 
increases  monotonically  from  a  value  of  zero  to  a  value  of  unity  in  the  time  it  takes  light 
to  traverse  the  diameter  of  the  charged  particle.  In  other  words,  7j(<)  is  close  to  a  unit  step 
function.  It  is  this  function  r}{s)  in  (1)  or  T]{t)  in  (2)  that  had  been  previously  overlooked 
and  that  Yaghjian  found  in  his  careful  derivation  of  (1)  and  (2).  This  unobtrusive  function 
completely  eliminates  the  noncausal  pre-acceleration  that  has  plagued  the  original  equations 
of  motion  ever  since  they  were  derived  by  Lorentz  and  Abraham  about  100  3'ears  ago. 

Although  (2)  is  a  rather  complicated  equation  of  motion,  our  preliminary  work  with 
particles  confined  to  rectilinear  motion  shows  that  closed-form  solutions  are  indeed  possible. 
In  particular,  for  the  velocity  u  in  the  x  direction,  that  is,  u  =  ux,  we  can  find  a  general 
solution  to  (2)  by  making  the  change  of  variables 

I 

sinh  (V/c)  =  c-yu  ,  dr  =  dZ/q  (3) 


with 

to  convert  (2)  to 


7  =  ^  =  cosh  (V/c) 

Fext{r)  =  mV'(T)  - 


(4) 

(5) 


Equation  (5)  is  a  first  order  linear  differential  equation  for  V'(r). 
given  by 

r  0. 


V'(r)  = 


Its  solution  for  all  r  is 


r  <  0 
r  >  0 


(6) 


where  the  external  force  is  applied  at  t  =  0  and  is  zero  for  r  <  0.  Integration  of  (6)  over  time 
produces  the  solution  for  the  velocity  of  the  charge  that  is  zero  for  r  <  0  and  continuous  for 
all  r,  even  across  r  =  0,  as  long  as  FeiK''')  continuous  or  has  a  finite  jump  across  r  =  0. 
In  other  words,  the  inclusion  of  the  ry  function  in  the  equation  of  motion  has  eliminated  the 
pre-acceleration  from  the  solution  to  the  original  equation  of  motion  without  introducing 
false  discontinuities  in  velocity  across  r  =  0  or  spurious  delta  functions  and  their  derivatives 
at  r  =  0. 


1.2  High-Frequency  Diffraction 

In  the  area  of  high-frequency  diffraction,  Dr.  Yaghjian  has  been  the  principal  investigator 
during  the  past  ten  years  responsible  for  the  development  of  incremental  length  diffraction  co¬ 
efficients  (ILDC’s)  that  have  been  incorporated  into  production  computer  codes  delivered  to 
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the  Air  Force.  This  highly  successful  work  began  with  his  discovery  that  three-dimensional 
(3-D)  ILDC’s  could  be  found  directly  from  the  two-dimensional  (2-D)  far  fields  of  planar 
canonical  scatterers  [6].  Shore  and  Yaghjian  immediately  applied  this  discovery  to  finding 
ILDC's  for  cracks  and  ridges  in  metal  surfaces,  and  this  led  to  the  first  accurate  prediction 
of  cross  polarized  fields  from  reflector  antennas  with  cracks  between  their  panels  [7],  Sub¬ 
sequently,  ILDC’s  were  found  for  truncated  wedges  [8],  and  both  the  truncated- wedge  and 
crack-ridge  ILDC’s  have  been  included  in  the  Air  Force  codes  for  computing  scattering  from 
large  bodies.  (The  truncated-wedge  work  won  an  IEEE  Best  Paper  Award  at  the  1996  IEEE 
International  Symposium  on  Antennas  and  Propagation.) 

Moreover,  Yaghjian  has  derived  a  theorem  that  removes  the  planar  restriction  on  2-D 
canonical  scatterers  used  to  find  3-D  ILDC’s.  With  this  restriction  removed,  he  has  been 
able  to  derived  (in  1996  and  1997)  the  first  shadow-boundary  ILDC’s  that  can  be  used  to 
significantly  enhance  the  accuracy  of  bistatic  as  well  as  monostatic  scattering  computed  from 
physical  optics  (PO)  codes.  These  shadow-boundary  ILDC  s  can  be  written  in  the  form 

dEssir',  e,  (»)  =  dB’ruir',  0,  d>)  +  9,  #)  +  0, «  +  «,  d>)  0) 


where 
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sin  uq 

The  TM  and  TE  shadow-boundary  electric  fields  required  on  the  right  hand  sides  of  (8)-(ll) 
are  given  in  [9].  Dr.  R.A.  Shore  of  Rome  Laboratory  tested  these  ILDC’s  by  applying  them 
to  compute  the  high-frequency  scattering  from  a  perfectly  conducting  sphere  of  radius  a 
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illurranated  by  a  plane  wave.  The  corrections  computed  from  the  incremental  fields  were 
added  to  the  PO  fields  of  the  sphere,  and  the  results  are  shown  in  Figures  1  and  2  for  the 
normalized  bistatic  scattering  cross  section  in  the  E-  and  H-planes  of  a  sphere  with 

ka  -  60.  Figures  1  and  2  demonstrate  that  the  incremental  fields  added  to  the  PO  fields 
produce  a  scattering  cross  section  much  closer  to  that  of  the  exact  eigenfunction  solution  of 
the  sphere  than  that  of  the  PO  solution  alone,  especially  for  large  bistatic  angles.  (The  exact 
and  PO+ILDC  fields  are  practically  indistinguishable  over  the  range  of  0  to  90  degrees  not 
shown  in  Figures  1  and  2.) 

The  combined  work  of  Yaghjian  and  Shore  in  determining  shadow-boundary  ILDC’s  was 
published  in  1997  as  an  invited  paper  to  a  special  issue  of  Radio  Science  on  electromagnetic 
theory  [9].  We  have  been  asked  by  Dr.  V.  Oliker  of  MATIS  to  help  implement  these  shadow- 
boundary  ILDC’s  in  the  near  future  into  the  production  high-frequency  computer  codes 
which  they  continue  to  develop  for  the  Air  Force. 


1.3  Surface  Integral  Equations 

In  the  area  of  surface  integral  equations,  Yaghjian  made  substantial  progress  in  1997  toward 
solving  the  problem  of  the  high  condition  numbers  one  encounters  when  applying  the  electric 
field  integral  equation  (EFIE)  to  perfect  conductors.  These  high  condition  numbers  increase 
the  number  of  required  iterations  in  an  iterative  solution,  and  thus  they  severely  limit  the 
size  and  accuracy  obtainable  using  the  EFIE.  To  quote  Professor  W.C.  Chew  from  his  March 
1997  review  article  in  the  IEEE  Transactions  on  Antennas  and  Propagation  (Fast  Solution 
Methods  in  Electromagnetics) 

Therefore,  preconditioning  techniques  for  reducing  the  required  number  of  iter¬ 
ations  in  iterative  methods  are  urgently  needed  in  solving  electromagnetic  wave 
scattering  problems.” 

Toward  this  end,  Yaghjian  was  able  to  determine  the  root  cause  of  the  EFIE  condition 
number  increasing  with  patch  density.  He  also  found  a  very  simple,  practical  preconditioner 
to  stabilize  the  condition  number  of  the  EFIE  applied  to  2-D  scatterers.  This  preconditioner 
which  was  presented  at  the  1997  lEEE/URSI  International  Symposium  in  Montreal  Canada’ 
is  presently  being  tested  m  the  2-D  EFIE  computer  codes  of  Boeing-McDonnell  Douglas  (St’ 
Louis,  MO)  and  Lockheed-Martin  (Palmdale,  CA).  »  v  • 

To  understand  why  and  hoW  this  EFIE  preconditioner  works,  it  helps  to  begin  with  the 
well-conditioned  magnetic-field  integral  equation.  Using  the  theorems  on  integral  equations 

from  [10],  It  was  proven  m  [11]  and  confirmed  numerically  that  the  spectral  condition  number- 
lor  the  magnetic-field  integral  equation  (MFIE) 

fi  X  fr,„e(r)  =  -hx  j>  K(r')  x  W(r,  r')dS'  +  (12) 

c  2 
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Figure  1:  Plot  of  PO,  PO+ILDC,  and  exact  normalized  bistatic  scattering  cross  section  of 
a  perfectly  conducting  sphere  (E-plane,  ka  =  60). 


Figure  2:  Plot  of  PO,  PO+ILDC,  and  exact  normalized  bistatic  scattering  cross  section  of 
a  perfectly  conducting  sphere  (H-pIane,  ka  =  60). 


(modified  by  a  ‘^dual  r  ’ 

value  as  the  number  of  patches  persXre^wTveWhT"^''"^  fUj)  approaches  a  finite 
conducting  scatterer  of  fixed  electrical  size  Cn  ^PP^aches  infinity  for  a  perfectly 

to  solve  the  dual-surface  MFIE  by  the  conjugTte'^era?’  f  ^  iterations  required 

scherne  also  approaches  a  constant  as  theTatch  size  f  "  similar  iterative 

tW  fllj,  u^n!  Dulse  h  f  '''  approaches 

that  the  number  of  conjugate  gradient  iterZn.  I  ^^tching, 

cu  e  of  surface  area  A  increases  as  (A/A^y/^  whe/eT^'^^tVi^^^  scattering  from  a 

time  required  to  iteratively  solve  the  MFIE  applied  to  a  /  CPU 

scatterer  can  be  approximated  by  the  propor^fonality  three-dimensional  (3-D) 


oc  {A/A'^yl'^Mi 


(13) 


where  TV  is  the  total  number  of  unknowns, 

for  the^rU  cu™i‘k  ''‘“'‘^e-field  integro-differential  equation  (EFIE)  (12] 

£;«(r)  =  X  I^K^ry)  -  ( v; .  K)VV(r,  („) 

has  a  spectral  condition  number  tbaf  ir,.- 
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method  and  a  strarghtforward  foment  nrethod  ^retS  ‘f*'" 

least  as  large  as  Discretization  scheme,  will  in  general  be  at 

Tboc  (f1/A2)'/2^AVA5)M2 
where  A5  is  the  uniform  patch  area  Pnr  =  f  u 

the  CPU  time  in  (IS)  needed  to  solve  the  EFEirTF'' wavelength 
t.me  .n  (13)  needed  to  iteratively  solve  tL  Mm  u  ^  *■"  '"S-  than  the  CPU 

reducing  the  iterative  solution  time  to  that  <r;,r  •  /,  ox  r  density  increases,  thereby 
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n  X 


,  \  -  -y  — ** 

X  [/  +  /  K,  •  vy'Hi^Uc' 
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and 


(17) 

where  H^\k\p  -  p'\)  is  the  Hankel  function,  K,(p)  is  the  current  tangential  to  the  curve 

Q  ZLd^nTe  >="«'*'  “>=  «lf-i-reme„t 

luZ  f"  conventional  TE  integro-differential  equation  correspond, ng  to 

14)  to  the  pure  integral  equation  in  (16)  (containing  no  derivatives  of  current  in  Ler 
to  allow  a  clearer  comparison  with  the  TE  MFIE  obtained  from  (12).  Of  course  the  TM 
electrm-field  Integra  equation  in  (1 7)  is  identical  to  the  conventional  TM  eleclric-field  i„tegr“ 
equa  ion.  Its  kernel  Ff„  (Ji|p  -  p'|)  ,s  fes  singular  than  the  kernel  of  the  TE  electric-field 
1  tegral  equation  (16).  and  it  can  be  solved  simply  and  straightforwardly  using,  forexlmpt 

TE  9^  matching.  We  have  also  solved  the  EFIE  in  (16)  for  the 

functions  and  point  matching.  Reasonable 
accuracy  was  found  for  the  current  using  as  few  as  5  increments  per  wavelength. 

(  4),  (16)  and  (17)  are  applied  to  closed  perfect  conductors,  they,  like  the  MFIE  inYl9) 
will  encounter  spurious  resonances.  For  such  closed  scatterers  ’ 

witf  orT;to™nnh?r:: . 

expressed  concisdTai’”’”  representation  of  (14),  (16)  or  (17) 

mi<]  =  I£,xJ  (18) 

V4  from  th.  diagona, 

about  A/4  of  the  diagonal  elements.  Then  (18)  is  transformed  to  the  set  of  equations 


|Z.-ZJ[A'I  =  (Z, 


(19) 


multiwavelength  bodies,  the  co^^puter  t,mX“^d  fo 
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to  a  small  fraction  of  the  total  computer  time  in  (13).  Moreover,  this  simple  banded-matrix 
preconditioning  converts  the  matrix  representation  (18)  of  the  electric-field  integral  equation 
(14),  (16)  or  (17)  into  a  matrix  representation  (19)  that  is  as  well  conditioned  as  the  MFIE, 
and  that  can  be  solved  iteratively  in  the  CPU  time  Tm  given  in  (13). 

These  conclusions  have  been  confirmed  by  solving  (18)  and  (19)  iteratively  for  TE  and 
TM  scattering  from  2-D  strips  with  widths  that  varied  from  one  to  forty  wavelengths.  Ta¬ 
bles  1  and  2  below  show  spectral  condition  numbers  and  numbers  of  iterations  for  the  TE 
scattering.  The  results  for  TM  scattering  are  similar  but  with  the  number  of  iterations  in 
the  preconditioned  cases  growing  more  slowly  than  the  (stripwidth/wavelength)^/^  growth 

exhibited  in  Table  2  for  the  TE  scattering. 
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Table  1.  TE  scattering  from  5- wavelength  strip 


Increments 

per 

Wavelength 

Number  of  Iterations 
without  Preconditioning 
(10“^  residual) 

Number  of  Iterations 
with  Preconditioning 
(10“®  residual) 

5 

22  (5.1) 

12  (2.7) 

10 

48  (11.1) 

11  (2.4) 

15 

72  (16.9) 

12  (2.7) 

20 

96  (22.7) 

13  (2.9) 

60 

did  not  converge  (57.4) 

13  (3.0) 

(Spectral  condition  numbers  given  in  parentheses) 


Table  2.  TE  scattering  from  strips  of  different  widths 
using  10  increments  per  wavelength 


Strip  Width 
in 

Wavelengths 

Number  of  Iterations 
without  Preconditioning 
(10“®  residual) 

Number  of  Iterations 
with  Preconditioning 
(10“®  residual) 

5 

48  (11.1) 

11  (2.4) 

10 

70  (15.3) 

14  (3.3) 

20 

101  (21.3) 

21  (4.6) 

40 

143  (29.8) 

30  (6.5) 

(Spectra!  condition  numbers  given  in  parentheses) 
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2  Principal  Publications  in  1997 

Books 

1.  T.B.  Hansen  and  A.D.  Yaghjian,  Plane-Wave  Theory  of  Time-Domain  Fields:  Appli¬ 
cation  to  Near-Field  Scanning,  New  York;  lEEE/Oxford  University  Press,  to  appear 
in  1998. 

Book  Chapters 

1.  A.D.  Yaghjian  and  T.B.  Hansen,  “Theorems  on  Time-Domain  Far  Fields,”  in  Ultra- 
Wideband,  Short-Pulse  Electromagnetics  3,  New  York:  Plenum  Press,  1997,  165-176 
(Invited). 

Journal  Articles 

1.  A.D.  Yaghjian,  R.A.  Shore  and  M.B.  Woodworth,  “Shadow-boundary  incremental 
length  diffraction  coefficients  for  perfectly  conducting  smooth,  convex  surfaces,”  Radio 
Science,  vol.  31,  1681-1695,  November- December  1996  (Invited). 
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radiation,”  Wave  Motion,  vol.  25,  1997. 

5.  A.D.  Yaghjian,  T.B.  Hansen  and  A.J.  Devaney,  “Minimum  source  region  for  a  given 
far-field  pattern,”  IEEE  Trans.  Antennas  and  Propagat.,  vol.  45,  911-913,  May  1997. 

6.  A.D.  Yaghjian,  “Banded-matrix  preconditioning  for  electric-field  integral  equations,” 
Digest  of  IEEE  Antennas  and  Propagat.  Soc.  Int’l  Symp.  (Montreal  Canada),  vol.  3, 
1806-1809,  July  1997. 

Citations  in  Science  Citation  Index 

Publications  authored  and  co-authored  by  Dr.  Yaghjian  have  received  an  average  of 
40  citations  per  year  for  the  past  10  years.  In  1996,  the  last  full  year  of  Science 
Citation  Index  published  to  date,  his  publications  received  46  citations.  For  every 
author  who  cites  a  publication,  there  are  usually  many  other  people  who  have  used 
that  publication  in  their  work  but  who  have  not  published  their  work.  Thus,  this  large 
number  of  citations  each  year  indicates  a  considerable  transfer  of  technology. 
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Appendix  A:  Derivation  of  Dispersion  Relations 


Causality  requires  that  the  time-domain  constitutive  parameters  are  zero  for  t  <  0.  W^e  can 
use  this  fact  to  relate  the  real  and  imaginary  parts  of  each  frequency-domain  constitutive 
parameter  by  means  of  the  Hilbert  transform.  Consider,  for  example,  the  frequency-domain 
complex  electric  susceptibility  xZ,  which  can  be  written  in  terms  of  the  Fourier  transform  of 
the  time-domain  electric  susceptibility 


Xu 


1  1  .  -t-co 

-  /  i  lim  !  i 

27r  /  27r  y  ^  27r  J  jtj 


(20) 


The  last  equality  in  (20)  can  be  proven  valid  for  <5  >  0  using  the  theorem  in  [18,  vol.2,  sec.225] 
that  allows  the  limit  to  be  brought  inside  the  integral  if  /+“  <  oo-  Assume  that 

Xl  is  a  Holder  continuous  function  of  cu,  and  for  sufficiently  large  |cuj,  ]x*  [  <  Cl\u\°‘,  where 
C  and  a  are  positive  constants  [19,  sec.43].  Then  one  can  apply  the  convolution  theorem  to 
the  last  Fourier  transform  in  (20),  use  the  result 


1  j. 

■-  - 

2w  J  |t| 


ILO 


(21) 


rewrite  the  convolution  integral  as  a  Cauchy  principal  value  integral,  and  take  the  limit  as 
^  0  (using  [18,  voL2,  sec.225]  and  the  fact  that  Xw  is  Holder  continuous),  to  obtain 


— oo 


(22) 


where  the  bar  through  the  integral  sign  in  (22)  denotes  the  Cauchy  principal  value  integration 

[20,  p.368].  Equating  the  real  and  imaginary  parts  of  (22),  one  obtains  the  “dispersion 
relations” 


-OO  0 


TT  J  i/  —  cu  7r 


“OO 


0 


(24) 


Equations  (23)-(24)  show  that  the  real  and  imaginary  parts  of  xl  are  Hilbert  transforms  of 
each  other  [20,  p.372].  In  fact,  these  dispersion  relations  could  be  obtained  directly  from 
the  theory  of  Hilbert  transforms,  given  that  the  first  equation  in  (20)  implies  that  xt  is 
an  analytic  function  in  the  upper  half  of  the  complex  w-plane.  Here,  however,  we  have 
derived  the  dispersion  relations  (23)- (24)  straightforwardly  without  relying  on  the  theory  of 
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Hilbert  transforms  or  integrating  in  the  complex  w-plane.  In  addition,  we  have  determined 
the  conditions  on  the  constitutive  parameters,  in  this  case  assure  the  existence  of 

the  Hilbert  transforms  and  the  validity  of  the  dispersion  relations. 

The  corresponding  dispersion  relations  hold,  of  course,  for  the  magnetic  susceptibility 
and  conductivity  provided  they  too  are  Holder  continuous  functions  of  w  that  approach 
zero  as  Cl\uj\'^  or  faster  (a  >  0)  as  |a;)  oo.  For  frequency-domain  functions  like  the 
complex  permittivity  =  6^  -f  ie"  =  e(l  -i-  x'),  which  approaches  the  free-space  value 
e  as  leu]  -4  oo,  the  dispersion  relations  (23)-(24)  can  be  applied  to  -  e)  to  yield  the 
“Kramers-Kronig  dispersion  relations”  [21,  p.311] 


—  OO  —CO 


“OO  — OO 


Formula  (23)  (or  (25))  is  especially  useful  because  approximate  values  of  the  odd  function 
Xw  ^  such  that  x^^^  ^  0  foi'  tu  >  0,  can  be  inserted  into  (23)  to  obtain  a  permissible  even 
function  xf  ■  However,  appi’oximate  values  of  the  even  function  x‘J  inserted  into  formula 
(24)  may  yield  an  odd  function  xV  does  not  satisfy  the  condition  xV  ^  0  for  u;  >  0 

[22,  p.281]. 

Each  of  the  dispersion  relations  (23)  and  (24)  (or  (25)  and  (26))  holds  exactly  only  if 
the  imaginary  and  real  parts,  respectively,  of  the  constitutive  parameter  are  known  and 
integrated  over  all  frequencies.  Approximate  dispersion  relations,  corresponding  to  (23)  and 
(24)  (or  (25)  and  (26)),  can  often  be  found,  nonetheless,  that  produce  accurate  values  of 
the  real  and  imaginary  parts  of  the  constitutive  parameter  from  the  local  variation  of  its 
imaginary  and  real  parts,  respectively  [23]. 
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Appendix  C:  Minimum  Source  Region  for  a  Given  Far-Field  Pat¬ 
tern 


A  number  of  authors,  including[24,  th.26]  and  [25,  sec.4.3],  have  proven  that  the  far-held 
pattern  produced  by  sources  in  a  volume  of  finite  extent  is  an  entire  analytic  function  of 
the  complex  spherical  angles  (0,  4>).  In  addition,  Mfiller  [24,  th.26,27,29]  has  proven  that  the 
sources  of  a  given  analytic  far-field  pattern  must  be  contained  in  a  sphere  with  a  radius  at 
least  as  large  as  R.  The  value  of  R  is  determined  from  the  spherical  expansion  coefficients 
of  the  far-field  pattern,  and  for  many  far-field  patterns  the  value  of  R  is  greater  than  zero, 
regardless  of  the  position  of  the  origin  chosen  for  the  spherical  expansion.  Outside  the  sphere 
of  radius  R,  the  field  of  the  given  analytic  far-field  pattern  is  also  an  analytic  function  of  the 
spherical  coordinates,  and  this  function  can  be  uniquely  continued  analytically  to  a  minimum 
source  region  V^in  located  inside  the  sphere  of  radius  R.  (This  can  be  done,  in  principle,  by 
expanding  the  field  in  a  series  of  Bessel  spherical  wave  functions  that  converges  uniformly  up 
to  the  surface  of  the  minimum  source  region  [24,  ch.3].)  Outside  K„,„,  which  can  be  singly 
or  multiply  connected,  the  field  obeys  the  homogeneous  wave  equation.  In  addition,  Knm  is 
the  smallest  possible  source  region,  in  that  the  sources  of  the  given  far-field  pattern  must 

extend  at  least  as  far  as  the  surface  of  or  else  the  field  could  be  further  analytically 
continued  inside  the  surface. 

Here,  we  derive  a  simple  method,  based  on  the  plane- wave  spectrum  formulation,  for 
determining  the  minimum  convex  closed  surface  in  which  the  sources  of  a  given  acoustic 
or  electromagnetic  far-field  pattern  can  be  located,  that  is,  the  minimum  convex  surface 
enclosing  the  minimum  source  region  Unlike  the  value  of  R  in  the  method  based 

on  a  spherical  wave  expansion,  the  shape  and  size  of  found  here  from  the  plane-wave 

spec  rum  is  independent  of  the  position  of  the  origin  of  the  coordinate  system  in  which  the 
plane-wave  spectrum  is  defined. 

We  begin  with  a  result  derived  in  our  book,  namely,  that  the  far-field  functions  FJK,  k^) 
for  sources  in  a  volume  U  of  finite  extent  satisfy  the  asymptotic  relation 


hi 


lim  =  0, 


z  >  z. 


for  all  ky/kx 


(27) 


where  z,  equals  the  maximum  (more  precisely,  the  supremum)  value  of  z  for  which  the 
source  function  ^(t,  y,  2)  is  nonzero.  Here,  can  be  a  negative,  zero,  or  positive  real 

number,  depending  on  the  position  of  the  origin  of  the  xyz  coordinate  system  with  respect 
to  the  source  region. 

An  infinite  number  of  diflPerent  source  distributions  can  radiate  the  same  far  field,  and 
thus  produce  the  same  field  everywhere  outside  the  source  region.  For  example,  any  number 
of  spherically  symmetric  source  distributions  produce  the  same  constant  far-field  pattern. 
Thus,  depending  on  the  particular  source  distribution  qUx,y,z),  (27)  may  hold  for  z  >  zi 
with  zi  <  z,.  Moreover,  unless  the  source  distribution  can  be  zero  everywhere  (and  thus 
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not  radiate  at  all),  there  must  be  a  minimum  for  which  the  limit  in  (27)  is  zero.  Letting 
^min  he  the  maximum  value  that  z  attains  in  so  that  z  ^  ^min  when  v  G  Vmin^ 

the  minimum  zi  is  equal  to  ;Zrnm-  To  prove  this,  assume  that  zj  were  less  than  so 
that  the  plane  z  z^  cuts  through  the  minimum  source  region  Vruin^  Then  the  spectrum 
ky)  =  (i/ j)Fu,{k^^ky)  would  decay  exponentially  according  to  (27)  for  z  >  zi,  and  the 
field  given  by  the  plane-wave  expansion 

J  +00  +00 

J  T^{kr,ky)e'(>‘^^+^yy+'^^Uk^dky  (28) 

—  c?o  —oo 

would  satisfy  the  homogeneous  wave  equation  for  z  >  z\  and  be  an  analytic  function  of  the 
spatial  coordinates  for  z  >  z\}  Consequently,  z  =■  z\  cannot  cut  through  the  minimum 
source  region,  and  the  minimum  value  of  zj  is  therefore  equal  to  Zm,„. 

Hence,  for  z  <  Zjnini  fhe  limit  in  (27)  is  indeed  not  equal  to  zero  at  least  for  one  value  of 
ky/kx  as  I7I  =  \kl  ky  —  —*■  00,  so  for  any  far-field  function  generated  by  sources  in  a 

volume  of  finite  extent,  we  have 


limsup|F,^(ifcr,jfcj,)je~l^l^  =  (  ^  ^  ^rnin,  iox&Wkylkx 

I7I— *-00  (  00,  z  <c  ^rn.in>  lor  at  Icast  one  kyjk^.. 


(29) 


The  infinite  value  occurs  on  the  lower  right  side  of  (29)  because  if  the  value  of  the  limit  is 
nonzero  for  all  z  <  Zmin,  the  exponential  factor  on  the  left  side  of  (29)  demands  that  its 
value  be  infinite  for  all  z  <  z,„.„.  The  “limit  superior”  [27,  p.68]  is  necessary  in  (29)  to  allow 
for  oscillating  values  of  F^{kx,  ky)  as  I7!  00. 

The  condition  (29)  can  be  rewritten  as  follows  in  terms  of  a  far-held  pattern  ^u,(^,  <p) 


limsup  ^  >  ^min,  for  all  4> 

fl— f-ioo  [00,  z  <  Zn{n,  for  at  least  one  ^  . 


(30) 


where  we  have  used  the  fact  that 


<^)  =  Fu,(k  cos  <j)  sin  6,  k  sin  (f)  sin  6)  (31) 

for  0  <  (^  <  27r;  0  <  ^  <  7r/2  and  6  =  7r/2  —  ia  with  0  <  a  <  00. 

Now,  suppose  we  are  given  an  analytic  far-field  pattern  and  have  determined  the  value 
of  Zmin  such  that  (30)  IS  satisfied.  The  minimum  source  region  for  this  far-field  pattern 
must  extend  in  the  z  direction  as  far  as  z  =  Zxn,in-  If  (30)  holds  for  all  orientations  of  the  z  axis 

^To  prove  that  the  field  in  (28)  satisfies  the  homogeneous  wave  equation  for  2  >  zi,  apply  the  operator 
to  the  right  side  of  (28)  and  interchange  the  integrations  and  differentiations  for  2  >  zj  as  allowed  by  the 
exponential  decay  of  the  integrand  [26,  sec.6.2,3].  Similarly,  it  follows  that  the  field  in  (28)  is  an  analytic 
lunction  of  X,  t/,  and  z  in  domains  that  include  the  real  x,  t/,  and  z  axes. 
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Figure  3;  Minimum  convex  closed  surface  formed  by  the  locus  of  points  for  all 
orientations  of  the  z  axis. 

(with  Zynin  and  depending,  in  general,  on  the  orientation),  then  the  corresponding 

spectrum  Tuj{k^,  ky)  for  all  orientations  of  the  z  axis  can  be  inserted  into  (28)  to  generate  a 
function  Pu{‘f’}  that  satisfies  the  homogeneous  wave  equation  outside  the  minimum  convex 
closed  surface  Smin  formed  by  all  the  Zmin',  as  shown  in  Figure  3.  In  other  words,  5mm 
will  be  the  minimum  convex  surface  enclosing  the  minimum  source  region.  The  plane-wave 
representation  (28)  also  gives  us  the  value  of  p^(r)  and  all  its  derivatives  on  Smin  In 
particular,  it  gives  the  Huygens’  sources  (p^(r)  and  {d ! dn)p^{T))  that  generate  the  fields 
outside  Smin- 

Consequently,  we  have  proven  that  a  far  field  satisfying  the  conditions  (29)  or  (30)  for 
every  orientation  of  the  z  axis  can  be  generated  by  sources  that  lie  within  or  on  the  minimum 
convex  closed  surface  Smin  determined  by  the  locus  of  values  of  Zmin  for  all  orientations  of 
the  z  axis.  Moreover,  this  far  field  cannot  be  generated  by  sources  that  extend  to  a  radius 
that  is  less  than  the  radius  Rmin  of  the  sphere  that  circumscribes  Smin  (see  Figure  3).  If 
the  origin  of  this  circumscribing  sphere  is  chosen  as  the  origin  of  the  coordinate  system  that 
Muller  [24]  uses  to  expand  the  fields  in  spherical  waves,  then  this  Rmin  equals  the  minimum 
R  in  [24,  th.26,27].  For  other  origins,  Mfiller’s  minimum  R  is  greater  than  Rmin-  Of  course, 
the  sources  that  generate  the  far-field  pattern  <j>)  are  not  unique  and  can  be  chosen  to 
extend  outside  Smin- 
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For  many  far-field  patterns,  Smin  will  simply  be  a  single  point  in  space  so  that  Rmin  =  0. 
This  means  that  the  pressure  associated  with  <f>)  can  be  produced  by  a  convergent 
sum  of  spherical  multipoles  located  at  a  single  point  in  space  If  Rmin  >  0,  a  multipole 

expansion  (located  at  the  center  of  the  sphere  with  radius  Rmin)  of  the  pressure  Pa-(»’)  having 
the  far-held  pattern  ^ will  converge  for  radii  r  >  Rmin-  Moreover,  at  least  for  some 
spherical  angles  {e,(f>),  it  will  diverge  when  r  <  R^i^.  If  the  infinite  sum  of  spherical 
inultipoles  is  truncated  to  a  finite  mode  number  N,  the  resulting  approximate  field  will 
converge,  of  course,  for  all  r  >  0  but  will  have  extremely  high  reactive  fields  for  r  <  N/k 
[28], [29]. 

For  electromagnetic  fields,  the  scalar  acoustic  far  field  is  merely  replaced  by  the  vector 
electromagnetic  far  field,  and  (29)  and  (30)  are  replaced  by 


limsup|F<^(ib^,jtj,)|e“l^l^  =  I  >  ^min,  ior  al]  ky/k^ 

|7|-too  (  CO,  z  <  Zmim  for  at  least  one  ky/k^ 


limsup  = 


0, 

oo, 


^  ^  ^min^  for  all  (j) 

■2  <  2mtn,  for  at  least  one  <j)  . 


(32) 

(33) 


Example 

We  shall  now  use  of  the  acoustic  formula  (30)  to  determine  the  minimum  convex  surface 
Smin  for  the  far-field  pattern  of  a  circular-disk  acoustic  radiator  of  radius  a.  This  far-field 
pattern  as  a  function  of  spherical  angles  (0,^)  measured  with  respect  to  a  rectangular  xyz 
coordinate  system,  with  z  axi^  normal  to  the  circular  disk  and  passing  through  the  center  of 
the  disk  at  2  =  0,  is  given  by 


T-  /n  i\  A  J\{ka  sin  ^)  „  „ 

<j>)  =  ’  ,  0<e<n  (34) 

where  Jy{x)  is  the  first  order  Bessel  function,  and  /4  is  a  constant  independent  of  6  and 
<6.  For  this  particular  orientation  of  the  ir  axis,  the  far-field  pattern  is  independent  of  the 
azimuthal  angle  <j>. 

Substituting  the  far-field  pattern  (34)  into  (30)  shows  that  Zmin  =  0  for  this  z  axis  normal 
to  the  center  of  the  disk.  For  the  z  axis  pointed  in  the  opposite  direction,  Zmin  s-lso  equals 
zero,  because  sin(7r  -  6)  =  sin  that  is,  the  far-field  pattern  in  (34)  is  identical  to  the 
right  and  left  of  the  circular-disk  radiator.  Thus,  we  have  confirmed  that  the  sources  of  the 
far-field  pattern  (34)  can  all  lie  in  the  plane  of  the  circular  disk  that  radiates  this  far-field 
pattern. 

To  find  the  minimum  distance  that  the  sources  must  extend  from  the  center  of  the 
disk  in  a  direction  in  the  plane  of  the  disk,  rotate  the  z  axis  90°  about  the  y  axis  and 
rewrite  the  far-field  pattern  in  terms  of  the  new  spherical  angles  (0^,  (j>')  measured  with 
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respect  to  this  rotated  z  axis,  which  now  lies  in  the  plane  of  the  disk.  Specifically,  we  have 
sin^  =  (1  -  sin^^'cos^  and  (34)  shows  that 

0  <  o  <  oo  (35) 

where  we  have  also  used  sin(7r/2  —  ia)  =  cosh  a.  Since 
(1  —  cosh^  a  cos^  ~  iz|  cos  (f)'\  cosh  a, 
the  asymptotic  formula 

J\[±ix) 

xsJl-Kx  ’  ^ 

can  be  invoked  to  get 


Tuj{T^j2  -  ia,(^')  =  A- 


Ji  ^A:ayl  —  cosh^acos^^ 
yjl  —  cosh^  a  cos^  <j}' 


limsup  \T^{Trl2  -  ia, 


A  cosh  a  I  cos  I  /czsinha 

1  •  C'  C 

=  limsup - 

or-*oo  \/2'!:ka  (cosh  a(  cos  ^'()3/2 


(38) 


when  cos4>'  ^  0.  For  cos^'  =  0,  we  have  T^{-k12  -  iaj')  =  AJr{ka).  Combining  these 
results  shows  that 


limsup  =  I  ^  ^ 

9'— f-ico  (_  OO,  z  <  a,  for  at  least  one  (j)' . 


(39) 


Comparing  (39)  with  (30),  we  see  that  z„„„  =  a  for  any  orientation  of  the  z  axis  in  the 
plane  of  the  disk.  Therefore,  the  minimum  convex  surface  Smin  containing  sources  that  can 
produce  the  far-field  pattern  of  the  circular-disk  radiator  simply  coincides  with  the  actual 
infinitesimally  thin  physical  disk  of  radius  a.  Of  course,  the  same  far-field  pattern  could  be 
excited  by  sources  that  extend  outside  this  disk.  However,  the  above  exercise  has  proven 
that  the  physical  disk  comprising  the  circular  acoustic  radiator  is  the  smallest  source  region 
that  can  exactly  produce  its  far-field  pattern. 

In  summary,  when  an  exact  closed  form  expression  for  the  acoustic  far  field  is  known,  (29), 

or  equivalently  (30),  determines  the  minimum  convex  closed  surface  in  which  the  sources  of 
the  far  field  can  reside. 
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